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We show that isotope-exchange reactions between ground-state alkali-metal, alkaline-earth-metal, 
and lanthanide heteronuclear dimers consisting of two isotopes of the same atom are exothermic 
with an energy change in the range of 1-8000 MHz thus resulting in cold or ultracold products. For 
these chemical reactions there are only one rovibrational and at most several hyperfine possible 
product states. The number and energetics of open and closed reactive channels can be controlled 
by the laser and magnetic fields. We suggest a laser-induced isotope- and state-selective Stark shift 
control to tune the exothermic isotope-exchange reactions to become endothermic thus providing 
the ground for testing models of the chemical reactivity. The present proposal opens the way for 
studying the state-to-state dynamics of ultracold chemical reactions beyond the universal limit with 
a meaningful control over the quantum states of both reactants and products. 
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Ultracold molecules have recently become a new pow¬ 
erful tool for investigating chemistry at its fundamen¬ 
tal quantum limit [1], The production of an ultra¬ 
cold high phase-space-density gas of polar molecules in 
their absolute rovibrational ground state [2-4] allowed for 
groundbreaking experiments on controlled chemical reac¬ 
tions [5-7]. An unprecedented control over the quantum 
states of reactants both by selecting their internal states 
and by tuning dipolar collisions with an external elec¬ 
tric field in a reduced dimensionality were achieved. The 
next expected milestone towards the state-to-state ultra¬ 
cold controlled chemistry will be the measurement and 
control of the product-state distributions [8] . 

The quantum statistics [5] , quantum tunneling and re¬ 
flection from a centrifugal barrier [9], as well as the im¬ 
pact of external fields [7] can modify the rates for ul¬ 
tracold chemical reactions by many orders of magnitude 
because at ultralow temperatures the reactants can be 
prepared in a single quantum state and even a weak per¬ 
turbation can be comparable to or larger than the colli¬ 
sion energy [8] . At the same time, the energy released in 
a chemical reaction between 1 E + -state polar alkali-metal 
molecules in the ground rovibrational level, if energeti¬ 
cally allowed, is at least 12 K for KRb, or more for other 
molecules [10] . This means that untrappable products in 
many states with relatively large kinetic energy are ex¬ 
pected and thus prospects for observing quantum effects 
and control in product distributions are limited [11]. 

The rate constants for ultracold reactive collisions of 
the KRb molecules were shown to be universal, that 
is, independent of the short-range dynamics when all 
wave function amplitude transmitted to a short range 
is lost [9]. To avoid reactive losses, which could ob¬ 
struct the interesting applications in quantum simula¬ 
tions of many-body physics [12], the creation of polar 
alkali-metal molecules in the rovibrational ground state 


stable against chemical reactions emerged as another im¬ 
portant research goal [3, 4], However, recently it was 
pointed by Bohn and collaborators [13-15] that even 
chemically stable molecules can form long-lived collision 
complexes which can lead to the delayed three-body loss 
mechanism. This prediction based on the statistical ap¬ 
proach together with other open questions on the dynam¬ 
ics of the reaction complex, non-universal reactivity, and 
potential control of state-to-state transitions still require 
experimental verification. For the systematic investiga¬ 
tion of ultracold chemical reactions, model systems with 
a small number of entrance and exit channels, and signif¬ 
icant control over the quantum states of both reactants 
and products are of great importance. 

Here we propose ultracold isotope-exchange reactions 
between heteronuclear dimers—a prototype of atom- 
exchange reactions—as a promising system for investigat¬ 
ing models of the controlled chemical reactivity. We show 
that the isotope-exchange reactions between all heteronu¬ 
clear alkali-metal, alkaline-earth-metal, and lanthanide 
dimers are exothermic with the reaction energy change 
as small as 1 MHz thus resulting in only one rovibra¬ 
tional and at most several hyperfine possible product 
states with the kinetic energy in the cold or ultracold 
regime. We demonstrate that these exothermic chem¬ 
ical reactions can be tuned to become endothermic by 
employing light near-resonant with a molecular bound- 
to-bound transition in the heteronuclear reactant while 
off-resonant for homonuclear products, that is, isotope- 
and state-selective AC Stark shift control. Application 
of the magnetic field can bring an additional control over 
the number and energetics of open and closed reactive 
channels. The present proposal thus paves the way for 
studying ultracold chemical reactions beyond the univer¬ 
sal limit with a meaningful control over the quantum 
states of both reactants and products. 
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The AC Stark shift control over the entrance chan¬ 
nels and reaction complex was recently proposed for a 
bimolecular reaction [16] and experimentally realized for 
unimolecular [17] and atomic [18] processes. 

Isotope-exchange reactions. — We consider chemical 
reactions between the ground-state heteronuclear dimers 
consisting of two isotopes of the same atom 

n A m A + n A m A n A 2 + m A 2 . (1) 


They will lead to collisional losses only when energetically 
allowed. Within the harmonic oscillator approximation 
used to describe the rovibrational ground states, the en¬ 
ergy change for the reaction (1) is given by: 


A E„ 


up (M m - M n ) 2 
128/r 2 


( 2 ) 


where uip is the harmonic frequency of the ground state, fi 
is the reduced mass of the m A n A molecule, and M m — M n 
is the mass difference of the isotopes m A and n A. It 
means that all isotope-exchange reactions of the type (1) 
are exothermic and all heteronuclear dimers are chemi¬ 
cally unstable at ultralow temperatures. The energy re¬ 
leased in the process is much smaller than the vibrational 
spacing, A E ex <C wo, therefore the products are possible 
only in their vibrational ground state. If the reaction 
change is smaller than the lowest rotational excitation 
energy, A E ex < 2 Bq (Bq is the rotational constant), then 
products are possible only in their rovibrational ground 
state. Equation (1) predicts that smaller reaction energy 
is expected for heavier dimers with smaller harmonic fre¬ 
quency and smaller mass difference of involved isotopes. 

To give an accurate prediction of the energy change 
for the isotope-exchange reactions we have performed 
full calculations for selected alkali-metal and alkali-earth- 
metal dimers by subtracting the rovibrational energies of 
the ground-state reactants from those of products. The 
rovibrational energy levels have been calculated from the 
available accurate experimental potential energy curves 
of the X 1 E + electronic ground state of the Li 2 [19], 
I <2 [20], Rb 2 [21], Sr 2 [22] molecules and of the a 3 E + state 
of the Li 2 [19], K 2 [20], Rb 2 [21] molecules. The hyper- 
fine structure and mass-dependent adiabatic and QED 
corrections to the Born-Oppenheimer interaction poten¬ 
tials are negligible [23]. The expression (2) reproduces 
the exact values within only a few percent error, there¬ 
fore the reaction energy changes for the Yb 2 , Dy 2 , Eu 2 
dimers, for which the experimental data are not avail¬ 
able, were obtained with this equation and theoretically 
predicted harmonic frequencies [24-26]. 

The energies released in the isotope-exchange reaction 
between dimers in the rovibrational ground level of the 
4 E + electronic state are presented in Table I and take val¬ 
ues between 0.87 MHz for 173 Yb 174 Yb up to 8104 MHz 
for 6 Li 7 Li. For 85 Rb 87 Rb and 87 Sr 88 Sr they are 29.1 MHz 
and 5.0 MHz, respectively. The formation of products 


TABLE I. The energy change A E ex for the reactions 
2 m A n A —> m A 2 + n A 2 of heteronucelear dimers in the rovi¬ 
brational ground state of the A X E + or a 3 E + electronic state 
and the estimated Rabi frequency H of the laser field needed 
to suppress reactivity ( AEac = §A7? ex ) assuming the detun¬ 
ing A nm = 10T e . Bo and too are the rotational constants 
and harmonic frequencies, respectively, for the rovibrational 
ground state of the rn A n A dimer. 


dimer 

AE ex (MHz) 

n (MHz) 

B 0 (MHz) 

wo (GHz) 

X 1 E + 

6 Li 7 Li 

-8104 

8466 

21700 

10920 

39 k 40 k 

-55.3 

207 

1660 

2746 

39 k 41 k 

-214 

447 

1640 

2729 

85 Rb 87 Rb 

-29.1 

148 

663 

1720 

87 Sr 88 Sr 

-4.98 

55 

527 

1216 

86 Sr 88 Sr 

-20.2 

112 

530 

1219 

84 Sr 88 Sr 

-83.2 

240 

536 

1226 

173 Yb 174 Yb 

-0.87 

22 

282 

746 

i68 Yb i76 Yb 

-56.3 

189 

284 

750 

163 Dy 164 Dy 

-0.90 

24 

287 

769 

160 Dy 164 Dy 

-14.5 

97 

290 

772 

151 Eu 153 Eu 

-3.7 

47 

276 

680 

o 3 E+ 

6 Li 7 Li 

-1510 

1843 

8650 

1992 

39 k 40 k 

-12.7 

96 

770 

625 

39 k 41 k 

-49.1 

194 

761 

622 

85 Rb 87 Rb 

-6.87 

72 

315 

401 


is energetically allowed only in the rovibrational ground 
state. The energy of 20 MHz corresponds to 1 mK, there¬ 
fore the products of all considered reactions have kinetic 
energy in the cold or ultracold regime. If the external 
magnetic field can be used to restrict the total electron 
spin of the collision complex to its maximum value [27], 
then the low-energy isotope-exchange reactions can be in¬ 
vestigated between alkali-metal dimers in the 3 E + state, 
for which the energy change is smaller by a factor of 4-5 
as compared to the values for the 1 E + state [cf. Table I]. 

Stark shift control with laser field. — Even small ex¬ 
ternal perturbation can drastically modify the dynamics 
of the reaction (1) with such a small energy change. Un¬ 
fortunately, the DC electric field control is not feasible for 
the ground-state homoatomic molecules. Here, we pro¬ 
pose to use the isotope- and state-selective Stark shift 
induced by the cw laser field [28] to control the energet¬ 
ics of the reaction (1). 

The AC Stark shift for an effective two-level system in 
the weak-field regime (f2 <C A) is given by: 


AEac 


nn 2 A 
4A 2 + r 2 ’ 


( 3 ) 


where Q. is the Rabi frequency of the laser, A is the de¬ 
tuning, and T e is the natural line width of the transi- 
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FIG. 1. (Color online) (a) Scheme of the state-selective AC 
Stark shift control. Different levels corresponds to different 
isotope mixtures (energies are not in scale), (b) The AC Stark 
shift for different laser detunings vs. the laser field Rabi fre¬ 
quency fi. (c) Zeeman levels of 8 'Rb 2 vs. the magnetic field. 


tion. The proposed control scheme is based on finding 
the laser wavelength which is near-resonant with some 
bound-to-bound transition [18] in the heteronuclear re¬ 
actant while off-resonant in the homonuclear products 
[cf. Fig. 1(a)]. If A nm < A mm ,A nn , where A mn is the 
detuning in the n A m A molecule, then the significant AC 
Stark shift is experienced only by the heteronuclear dimer 
and modifies the energy change for the isotope-exchange 
reaction (1). Potential heating and losses can be avoided 
by using the detuning much larger than the transition 
line width, A nm F e . To employ the above control 
scheme in a realistic molecular system the vibrational 
spacing for the ground and excited states must be larger 
than the detuning, — E™ m A nm , and the excita¬ 
tion energy for different isotopomers must differ by more 
than the detuning, A nm . The lat¬ 

ter is guaranteed when the ground level is the vibrational 
ground state, whereas the excited level is the vibrational 
state from the middle of the spectrum, v' v" = 0. All 
these conditions can be met for the investigated dimers 
and only an accidental coincidence can result in the same 
Stark shifts for the hetero- and homonuclear dimers. 

To ensure sufficiently large Rabi frequencies, we choose 
a transition to an excited state rovibrational level for 
which the atomic transition at the dissociation threshold 
is dipole allowed and strong [29]. For alkali-metal atoms 
this is the 2 Si /2 — 2 P3/2 transition at wavelengths from 
670 nm for Li to 780 nm for Rb with the line width of 
36-38 MHz. For Sr and Yb, the 1 5o — 1 Pi transition at 
wavelengths 461 nm and 399 nm with the line widths of 
30 MHz and 29 MHz, respectively. For Dy and Eu, the 
5 /s~ 5 ( 8 , 1)9 and s S’ 7 / 2 — s Pg /2 transitions at wavelengths 
412 nm and 466 nm with the line widths of 32 MHz and 
30 MHz, respectively. 


The molecular excitation energies and line widths will 
be different but of the same order of magnitude as the 
above listed atomic ones. Therefore, to estimate the Rabi 
frequency H of the laser needed to turn the exothermic 
isotope-exchange reactions into the endothermic ones, we 
compare the AC Stark shift for the heteronuclear dimers 
with the reaction energy change, 2A E^t(fl) = AE ex , 
assuming the detuning A" m = 10r e . The obtained 
Rabi frequencies are collected in Table I and are between 
22 MHz for 173 Yb 174 Yb up to 447 MHz for 39 K 41 K and 
8466 MHz for 6 Li 7 Li. The last one is quite large but 
for all other dimers the proposed AC Stark shift control 
should be experimentally feasible. 

To confirm the above predictions based on the simpli¬ 
fied two-level calculations, we performed the full calcula¬ 
tions of the AC Stark shift control by diagonalizing the 
molecular Hamiltonians for the 85 Rb 87 Rb and 87 Sr 88 Sr 
dimers for which accurate electronic interaction poten¬ 
tials and transition dipole moments are known [21, 22, 
30, 31]. For 85 Rb 87 Rb ( 87 Sr 88 Sr) we chose the laser 
wavelength A l = 619 nm (441 nm) and the detuning 
A nm = 370 MHz (305 MHz) from the level with the line 
width r e = 37MHz (30.5 MHz). The Rabi frequency 
H = 160 MHz (62 MHz) is needed to suppress the reactiv¬ 
ity. The calculations with the full Hamiltonians confirm 
the simplified predictions. 

Control with magnetic field .— The 4 S + -state 
molecules have a zero electronic spin therefore the 
hyperfine splitting of the order of 1-100 kHz due to 
the nuclear spin-spin coupling occurs only in the case 
of two nuclei with a non-zero spin [23]. If at least one 
nucleus in the heteronuclear reactant has a non-zero 
spin then products in a few hyperfine states are possible. 
Application of the magnetic field lifts the degeneracy of 
the hyperfine states and can be used to control energies 
of exit channels associated with different Zeeman levels. 
For example, the Zeeman shift for 8 ‘Rb 2 in the magnetic 
field of 1000Gauss is up to 4.2MHz [cf. Fig. 1(c)]. 

Applications .— The proposed cold isotope-exchange 
reactions together with the control schemes provide a 
platform for testing models of the chemical reactivity. 

The diversity of stable isotopes for Sr, Yb, or Dy can be 
used to tune the energy change, Eq. (2). The atoms with 
several bosonic and fermionic isotopes give the unique 
opportunity for realizing isotope-exchange reaction be¬ 
tween both fermionic and bosonic dimers and observing 
the effect of the quantum statistics. 

If the reaction proceeds via the formation of an inter¬ 
mediate complex with a sufficiently long lifetime, whose 
dynamics renders the reactant and product channels sta¬ 
tistically independent, then the probability for the state- 
to-state transition P a -+p (related to reaction rate) can be 
described within the statistical quantum model [32] by: 

Pa—t/3 ~ PaPfi/^^Pj ) 

7 


( 4 ) 
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where pi is the capture probability of the complex for¬ 
mation for the collision in the channel i. In the limit of 
many more exit than entrance channels, the total rates 
for chemical reactions are predicted to depend solely on 
the capture dynamics in the reactant channels [11] and 
the universal limit of reactivity can be expected [9]. 

In the investigated systems, the rotational excitation 
can be used to control the number of exit rovibrational 
channels in the range from one to several by Raman exci¬ 
tation [33] or to several tens if excited by using an optical 
centrifuge [34], This can allow to investigate the transi¬ 
tion from the non-universal limit of the chemical reac¬ 
tivity to the universal one [35] in a single experimental 
setup. The vibrational excitation can bring the system 
directly into the regime of the large number of exit chan¬ 
nels and thus the universal limit of chemical reactivity [9] . 

The reaction energy change can continuously be tuned 
from tens of MHz to zero with the proposed AC Stark 
shift control scheme. This should allow to observe the 
impact of the quantum statistics, Wigner threshold laws, 
and quantum tunneling and reflection from the centrifu¬ 
gal barrier on the capture probabilities in the product 
channels and thus the quantum control of reaction rates 
by controlling the product states. The selective activa¬ 
tion/deactivation of higher partial waves and tuning the 
position of shape resonances in the exit channels across 
the entrance threshold should also have a dramatic effect 
on the reaction rates and product distributions. The es¬ 
timated heights of the centrifugal barrier for the p-wave, 
d- wave, and /-wave collisions between Sr 2 or Rb 2 dimers 
are 0.35 MHz, 1.8 MHz, and 5.2 MHz, respecitively. 

The comparison of the collisional losses when the chem¬ 
ical reactions are allowed with those when they are sup¬ 
pressed by the proposed control scheme should allow to 
verify the predictions based on the statistical models 
about the long-lived collision complexes formation and 
subsequent losses [13-15] in a single experimental setup. 
More specifically, it was predicted that the lifetime of the 
collision complex r is proportional to the density of the 
available rovibrational states of the complex p and in¬ 
versely proportional to the number of open channels N „, 
r ~ p/N 0 [14]. In the proposed experiments, the number 
of open channels can be controlled by rotational excita¬ 
tion or by using the combined Stark shift and magnetic 
field control with the number of open Zeeman levels from 
one to several. 

The investigated reactions proceed barrierlessly via the 
four-body complex stabilized by non-additive forces [10, 
27] . The complex dynamics is still an open question and 
proposed experiments should also verify in a controlled 
way the statistical independence of the reactant and 
product channels and the existence of spectator degrees 
of freedom not involved during the reaction, that is, the 
conservation laws of different angular momenta present 
in reacting dimers. The control of the isotope-exchange 
reactions between the spin-polarized triplet-state dimers 


in a magnetic field can potentially give some insight into 
the dynamics of the transitions to lower spin states me¬ 
diated by nonadiabatic spin-dependent couplings [36]. 

Experimental accesibility .— The preparation of ul¬ 
tracold dense gases of heteronuclear alkali-metal dimers 
in the temperature range of 1-100 pK should not be 
more challenging than the formation of heteronuclear 
40 K 87 Rb [2] and 87 Rb 133 Cs [3, 4] molecules which were 
produced in their rovibration and hyperfine ground state 
by using the magnetoassociation followed by the efficient 
stimulated Raman adiabatic passage (STIRAP). Many 
other ultracold heteronuclear Bose-Bose, Fermi-Fermi, 
and Bose-Fermi mixtures of alkali-metal atoms were real¬ 
ized and converted into Fesbhach molecules [37], includ¬ 
ing 85 Rb 87 Rb [38]. Actually, the formation of heteronu¬ 
clear alkali-metal dimers can potentially be facilitated by 
similar trapping conditions and laser wavelengths needed 
for different isotopes as compared to different atoms. 
The magnetoassociation of heteronuclear dimers from the 
highly magnetic lanthanide atoms has to face the problem 
of the high density of the interaction-anisotropy-induced 
Feshbach resonances [26, 39] but then the STIRAP sta¬ 
bilization should be feasible. The formation of heteronu¬ 
clear molecules out of closed-shell atoms is a more chal¬ 
lenging problem [40]. Nevertheless, the photoassocitive 
formation [41] of the homonuclear strontium [42] and yt¬ 
terbium [43], and heteronuclear ytterbium [44] dimers 
into weakly bound states were already demonstrated and 
a stabilization to the ground state can be expected. 

The total reaction rates can be measured by monitor¬ 
ing the time evolution of the average number density of 
trapped heteronuclear reactants [5-7]. The state-selected 
product distribution measurement should also be possi¬ 
ble for trappable ultracold homonuclear products [45]. 
The estimated heating rate due to the photon scattering 
at a near-resonant transition used in the AC Stark shift 
control of the investigated Rb 2 and Sr 2 dimers is of the 
order of 1 mK/s for A nm = 10r e and can further be re¬ 
duced by one to two orders of magnitude by increasing 
the used detuning A nm [46]. This should be sufficiently 
low to allow for milliseconds measuraments. 

Summary. — We have shown that the ultracold 
isotope-exchange reactions between heteronuclear dimers 
in external fields provide the ground for testing models of 
the chemical reactivity where for the first time a mean¬ 
ingful control over the quantum states of not only reac¬ 
tants but also products can be realized. These exother¬ 
mic reactions can be tuned to become endothermic by 
employing the laser-induced isotope- and state-selective 
Stark shift control and the magnetic field brings an ad¬ 
ditional control over the number and energetics of open 
and closed reactive channels. The demonstrated control 
schemes are also applicable to other isotope-exchange re¬ 
actions e.g. in ultracold mixtures of molecules and atoms. 
The proposed systems and required field strengths are all 
within the current experimental capabilities. 
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